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Abstract 

We calculate open heavy-flavor (HF) transport in relativistic heavy-ion collisions by applying a strong-coupling 
treatment in both macro- and microscopic dynamics (hydrodynamics and non-perturbative diffusion interactions). The 
hydrodynamic medium evolution is quantitatively constrained by bulk and multi-strange hadron spectra and elliptic 
flow. The heavy quark (HQ) transport coefficient is evaluated from a non-perturbative T-matrix approach in the 
Quark-Gluon Plasma (QGP) which, close to the critical temperature, leads to resonance formation and feeds into the 
recombination of heavy quarks on a hydrodynamic hypersurface. In the hadronic phase, the diffusion of HF mesons is 
obtained from effective hadronic theory. We compute observables at RHIC and LHC for non-photonic electrons and 
HF mesons, respectively. 
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1. Introduction 

A deconfined state of nuclear matter has been predicted by large-scale numerical simulations of Quantum Chro- 
modynamics (QCD) on the lattice HJ|2). Utilizing ultra-relativistic heavy-ion collision (URHIC) experiments, this 
new phase has been identified as a strongly coupled Quark-Gluon Plasma (sQGP) EH). With the advent of LHC, the 
endeavor of characterizing the properties of sQGP has entered a new stage |5]|6]. 

Heavy quarks, due to their large masses (m c =s 1.3 GeV for charm and - 4.2 Ge V for bottom), are produced in 
primordial hard collisions [7 1, and their number is expected to be conserved in the subsequent evolution of the medium. 
Through interactions with medium constituents, the spectra of HQs are modified, yet they may not fully thermalize. 
This makes them valuable probes of the hot and dense matter created in URHICs (8). The suppression of non-photonic 
electrons (e ± ) from HF decays measured in Au+Au( y/s^ = 200 GeV) at RHIC ll9l ITOl fTTTl is surprisingly large, and 
their elliptic flow iflOl [Tl I indicates a substantial collectivity of charm quarks. It is difficult to explain these results 
within perturbative QCD approaches Ifl2l [131 [141 [131 for HQ scattering in the QGP. They rather indicate the need for 
a non-pertubative treatment of HQ diffusion in the low- and intermediate-momentum region fl6l . 

Here, we report on applications of our recently developed non-perturbative framework for open HF transport in 
medium ifTTl to D- and B-meson observables at LHC, and predict upcoming RHIC data for flavor-separated e ± spectra. 

2. Theoretical Framework: HQ Transport and Hadronization 

The kinetics of HQs in the QGP can be described by Fokker-Planck dynamics, which in practice is simulated 
through relativistic Langevin diffusion (S). The required transport coefficient is the HQ thermal relaxation rate, 
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Figure 1: (Color online) Separate charm- and bottom-electron Raa (left panel) and V2 (right panel) in Au+Au 
(VsNN=200GeV). 

A(p, T). We employ the results from a thermodynamic r -matrix approach 1 18 1 for HQ-light quark scattering with 
input potentials taken as the HQ internal energy computed in thermal lattice QCD. While HQ-gluon scattering was 
treated perturbatively in Refs. lfT71[T9ll . we here adopt a recent update also computed in the 7"-matrix formalism [20|. 
The T-matrices exhibit resonant states close to threshold at temperatures T - 1 — 1.5 T c , which accelerates the HQ 
thermal relaxation by a factor of ~3-5 for the HQ-light quark scattering (factor ~2 for HQ-gluon) over leading order 
pQCD results [18|. After hadronization, the open HF mesons continue to diffuse in the hadronic fluid. The ther- 
mal relaxation rate of D mesons has been calculated by using elastic scattering amplitudes from effective hadronic 
theory EH . As summarized in Ref. fl9l . the charm diffusion coefficient D s = T /[m c oMp - 0,T)] is essentially 
continuous across the phase transition region, possibly indicative of "quark-hadron duality". The B-meson thermal 
relaxation rate used is obtained from that of D mesons by scaling the relevant mass, assuming they have the same 
spatial diffusion coefficient (in the large mass limit). 

The Langevin simulations of HQ diffusion in QGP are implemented into the space-time evolution of a background 
medium modeled by boost-invariant ideal hydrodynamics [22 1. The equation of state (EoS) was constructed from a 
parameterization of recent lattice-QCD data in the high temperature phase 1231 121. combined with a partial chemical- 
equilibrium EoS for the hadronic phase. A compact initial-density profile with pre-equilibrium flow was introduced 
resulting in fair fits to multistrange-hadron observables close to the transition temperature, T c = 170 Me V, and to 
bulk-hadron observables at 7\j n = 110 MeV. Further details of the implementation of the Langevin simulation can be 
found in Ref. ifTTl . 

At T c = 170 MeV, the heavy quarks are hadronized into HF mesons. We accomplish this by first applying reso- 
nance recombination E4l[T71l with thermal light quarks on the pertinent hydro-hypersurface. Remaining heavy quarks 
are treated with (5-function fragmentation (as used in the fits to baseline HQ spectra in p+p collisions). For a reliable 
recombination dynamics at low and intermediate transverse momenta, pj, it is crucial that the correct thermal equi- 
librium limit is respected. We have verified this in the presence of the full space-momentum correlations generated 
from the HQ Langevin dynamics and hydrodynamic flow of the thermal light quarks. The relative partition between 
recombination and fragmentation is determined by a coalescence probability which is obtained from the HQ T-matrix 
scattering rate, and thus consistent with the underlying transport coefficient [17]. At T c , the in-medium quark and 
meson masses are used: m q = 0.3 GeV, m c = 1.7 GeV, nif, = 5.2 GeV, mo = 2.1 GeV and mg = 5.6 GeV, with D- and 
B-meson widths of Fd,b =0.1 GeV. After hadronization, the Langevin simulation resumes for D and B mesons in the 
hadronic phase until hydrodynamic kinetic freezeout at 7\j n = 1 10 MeV. 



3. Open Heavy-Flavor Observables 

3.1. Flavor-separated non-photonic electrons at RHIC 

The measured e ± nuclear suppression factor, R\a, and elliptic flow, V2, for Au+Au ( ysNN=200GeV) collisions 
thus far are superpositions of open-charm and -bottom contributions ||9j [TO] [TT] . In anticipation of experimental 
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Figure 2: (Color online) D-meson Baa (left panel) and V2 (right panel) in Pb+Pb ( yfs^-2.76 TeV) collisions, com- 
pared to ALICE data [28, 29 1. Our results for charm quarks, D mesons after hadronization at T c and after hadronic 
diffusion at 7\i n , are shown separately. The bands indicate the prevalent uncertainties in each case: shadowing (64.0%- 
77.6%) for Baa and charm-quark coalescence probability (50%-90%) for v%. 




Figure 3: (Color online) B-meson Baa (left panel) and V2 (right panel) for Pb+Pb ( -^snn=2.76 TeV) collisions. No 
shadowing is assumed. The average pj of the CMS datum 1 30 1 on non-prompt J/tf/ has been rescaled by m^/ray/^. 

results for charm- and bottom-separated e ± [25 1, we present in Fig. [TJour pertinent predictions. The charm-e* Baa 
exhibits a flow bump at p e t 1 .5 GeV, which is a combined consequence of charm-quark thermalization in the QGP 
diffusion processes and the subsequent coalescence with thermal light quarks at hadronization 1171 . Note that, from 
a microscopic point of view, both processes are of identical origin (heavy-light quark r-matrix). At p e T 5 GeV, the 
charm e ± suppression reaches 0.2-0.3. In accord with approaching thermalization in the suppression factor, the charm- 
e* also acquires substantial elliptic flow, reaching ~8% at p e t =s 2 GeV. Bottom electrons suffer less suppression, at 
least below p e t — 5 GeV, due to the larger b-quark mass. We also observe the bottom-e* v'2 to saturate at a significantly 
larger p e t than for charm, as a result of the larger mb/m q ratio that is operative in picking up light-quark V2 through 
coalescence. 



3.2. D and B mesons at LHC 

Next, we apply our non-perturbative approach to the LHC. Toward this end, we employ FONLL pQCD ll26ll Do- 
meson and B-meson spectra in V* = 2.76 TeV p+p collisions as our baseline, which yield good agreement with the 
ALICE D spectrum Il27l . Assuming ^-function fragmentation and folding with EPS09 shadowing [28 1, we obtain 
our initial HQ spectra for the Langevin simulations. The hydrodynamic medium evolution model has been tuned to 
available hadron data at the respective centralities, reproducing fairly well charged-hadron multiplicities, spectra and 
V2 at 7kj n , and Q~ observables at T c . 

In Fig.|2]we display our calculated D-meson Baa and \>2 for central and semi-central collisions, respectively. The 
flow bump in the c-quark B A a at low p T is amplified via coalescence at the D-meson level. Diffusion of D mesons in 
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the hadronic phase slightly reduces the D-meson Baa at high pj, and increases the V2 by up to 20%. Our D-meson 
Raa appears to be systematically slightly too high but it reproduces the /j^-shape of the data [28] rather well, largely 
induced by the falling momentum-dependence of the non-perturbative c-quark relaxation rate [18. 20]. The deviation 
from the data possibly indicates a missing radiative energy-loss contribution. For the D-meson we note that c-quark 
diffusion through the QGP only accounts for about half of its final value. Coalescence with light quarks and D-meson 
diffusion in hadronic phase cannot be neglected. Again, the apparent underprediction of the ALICE data, especially 
at higher pj where the the momentum-dependence of the elastic c-quark relaxation rate drops significantly, might 
call for an inclusion of radiative contributions. We note that such an assessment is only possible in the absence of 
phenomenological K-f actors. 

The B-meson results, shown in Fig. [3] show a substantial suppression of R AA 0.4 at p T > lOGeV, while the Vz, 
reaching up to 5%, indicates less collectivity (and thus less thermalization) than for D's. The B-meson modifications 
also impact J/if/'s via a so-called "non-prompt" feeddown. In this context it is interesting to quote the suppression of 
non-prompt J/if/of the CMS collaboration 11301 . which is ~0. 37 ±0.08 at an average p 3 j^ 9.3 GeV, in good agreement 
with our B-meson result at a correspondingly somewhat larger parent momentum. 

4. Conclusions 

We have evaluated heavy-flavor probes at RHIC and LHC within a non-perturbative framework for bulk evolution 
and microscopic interactions in the diffusion process 1 17 1. Initial comparisons to e ± , D- and B-meson observables at 
RHIC and LHC are encouraging. 
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